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We determine the mass of the lightest qq scalar meson with 1=1 using the simulation with two 
dynamical Domain Wall Fermions. The conventional exponential fit of the scalar correlator is 
justified in this case giving the mass 1.58 ±0.34 GeV. In general the scalar correlator receives 
also the bubble contribution, which is the intermediate state with two pseudoscalar mesons. This 
contribution is sizable at light quark masses and has to be incorporated in the fit of the scalar 
correlator in order to extract the scalar meson mass. We provide predictions for the bubble con- 
tribution in Partially Quenched ChPT, Staggered ChPT and ChPT for Mixed quark actions. We 
find that the bubble contribution is significantly affected by the unphysical approximations that 
are employed in simulations. It can render the negative sign and unphysical effective mass 2M K 
in the scalar correlator with 7=1. 
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1. Introduction 

The Nature of the lightest observed scalar resonances is not revealed yet. There are two exper- 
imentally well established scalar resonances ao(980) and ao(1450) with isospin 1=1 below 2 GeV. 
It is still not clear which one of the two is the lightest du scalar state. This raises a further question 
whether «o(980) is perhaps a tetraquark if ao(1450) turns out to be the lightest du scalar state. 

The first issue could be settled with a determination of the lightest du scalar mass on the lattice. 
For this purpose the lattice simulations evaluate the scalar correlator 

C{t) = '£(0\d(x,t)u(x,t) w(0,0)d(0,0)|0) =Ae- mMt +B(t) + -- - . (1.1) 



If ao is the lightest state with 1=1 and J p = + , the correlator (LI) drops as e~ mMt at large t and 
determination of m a o is straight forward. Multi-hadron states with J p = + and 1=1 also propagate 
between the source and the sink and they often shadow the interesting part e~ m,at in the correlator 



(1.1). The most important multi-hadron state is the intermediate state with two pseudoscalars, we 
call it the bubble contribution B{t) and display it in Fig. [[]. In proper three flavor QCD the two- 
pseudoscalar states are 7rq,KK, 717]', while Tin is not allowed by Bose symmetry and conservation 
of J p and I G . In two flavor QCD the only state icq' is relatively heavy and therefore not so disturb- 
ing. The scalar correlator ( pTT| ) is dominated by B{t) in the simulations with light quark masses if 
Mp\ +Mp2 < m a Q which opens the decay channel ao — > P\P2- So the bubble contribution has to be 
incorporated in the fit of the scalar correlator ( |1.1[ ) in order to extract m a Q. 

In addition to the physical intermediate states P\P2, the bubble contribution incorporates also 
the effects of the unphysical approximations employed in the simulations. This allows the extrac- 
tion of m a Q via ( |lTT| ) even in this case. The unphysical effects were first attributed to the bubble 
contribution in case of quenched QCD [gp, where it describes the magnitude and the negative sign 
of the lattice correlator well [||, ||]. 
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Figure 1: The bubble contribution to the scalar correlator. Here Pi and Pi denote pseudoscalar mesons in 
the relevant version of ChPT. 

Here we present the analytical predictions for B(t) in partially quenched simulations [Q] and 
in simulations with staggered fermions or mixed-quark actions [gj]. The results are obtained at the 
lowest order in the appropriate version of Chiral Perturbation Theory (ChPT) and apply for the 
point-point correlators. 

We demonstrate the method of determining m a o in a simulation with two dynamical Domain 
Wall Fermions (DWF) [Q]. The resulting scalar meson masses from the dynamical and partially 
quenched correlators agree. 

The simulation with chiral (DWF) sea and valence quarks is presented in Section 2. Section 3 
considers the simulations with staggered sea and valence quarks. The simulations with staggered 
sea and chiral valence are considered in Section 4. We conclude in Section 5. 
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2. Simulation with two dynamical Domain Wall quarks 



Dynamical correlator with m va \ 

The simulation with two Domain Wall sea quarks and Domain Wall valence quarks [||] does 
not suffer from unphysical effects in contrast to the simulations that will be considered later. The 
conventional exponential fit C(t) = Ae~'" a0t ( |0| ) is justified in two-flavor QCD since the only 
intermediate state nrj' is relatively heavy and B(t) is small [ffl]. The resulting m a o for three different 
quark masses is presented in Fig. 0a, while the linear extrapolation to the chiral limit gives 



in 



dyn 
aO 



1.58 ±0.34 GeV . 



(2.1) 



This mass has large error-bar, but it is definitely above 1 GeV and favors ao(1450) as the lightest 
du scalar meson. 

(a) Dynamical (b) Partially Quenched 
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Figure 2: Mass of ao from a simulation with two dynamical Domain Wall quarks, a ~ 1.7 GeV and 
V = 16 3 x 32: (a) m a Q from dynamical scalar correlator with am va \ = am sea using exponential fit; (b) m a o 
from partially quenched scalar correlator with am va i J=am sea = 0.02 using fit (1.1) with B PQChPT (2.2). 



Partially Quenched correlator with m va i ^ 

The partially quenched scalar correlator in Fig. ^ is positive for m va i > m sea , whereas it is 
negative for m va i < m sea H]. The striking effect of partial quenching is attributed to the bubble 
diagram in Fig. |T} In order to compute the bubble diagram we need the coupling of the point scalar 
current to two pseudoscalars and the pseudoscalar propagators from Partially Quenched ChPT. The 
coupling is obtained from du = —d^f chPT /d^du and is equal to Bo = Af^/ (2m q ) in any version 
of ChPT at the lowest order. The bubble contribution for mo — > °° is B(t) = F.T. [B(p)]p=o with [||] 



B p Q ChPT (p)=2Bl^< 



1 



1 



1 



k 2 +M 2 



(2.2) 

It does not contain any free parameters, since pseudoscalar masses M and Bo = M 2 /(2m q ) are 
determined from the lattice data. The sum over the loop momenta k is performed over the allowed 
discrete momenta on the lattice. Fig. || shows that the bubble contribution is positive for m va i > 
m sea and negative for m va i < m sea like the lattice data. The negativity of the scalar correlator for 
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fHval < msea and large t is described well by the bubble contribution. Having understood the source 
of the unphysical effect of partial quenching, we extract m a o by fitting the partially quenched scalar 
correlator to ( |1.1| ) and display it in Fig. ||a. The linear extrapolation to the chiral limit gives 



rrf a Q = 1.51 ±0.19 GeV in agreement with the dynamical result (|2JJ) and with smaller error. 
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0.001 



-0.001 



-0.002 
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m.=0.01 , m=0.02 




Figure 3: The symbols present the partially quenched point-point lattice correlators from a simulation with 
two dynamical Domain Wall quarks. The lines present the corresponding B p ® chPT (t) (2.2). 



3. Simulations with staggered sea and valence quarks 



The scalar correlator was simulated using 2+1 staggered sea quarks and staggered valence 
quarks by MILC [§] and UKQCD [f|. The authors of [|, [§] were surprised to find the effective 
mass significantly below M n +M,j, although %r\ is the lightest state with / = 1 for light u/d quarks 
in proper QCD. 

This can be attributed to the taste breaking effects which enter the scalar correlator via the 
bubble contribution [g] in the Staggered ChPT of [g] 



B SChPT {p)=B^ -4 



(k 2 +M 2 Si ) 



[k+pf+Ml 3 (*2 +M 2 / )(F + lM2 f + !Mf,) 



3 lrl U, 3 lrl S,< 



+ 



1 16 

-Y 



i 



{k + p) 2 +M 2 v 
1 



a 2 5y 



(k 2 +Mj v ) 



(k 2 + M 2 v )(k 2 +M 2 v )(k 2 +M 2 , 
1 1 



(3.1) 

+ (V^A) 



+ ■ 



1 



{k + p) 2 +M 2 k 2 +M 2 {k + p) 2 +M 2 Sb k 2 +M 2 



with the notation from [Rp- The forth-root trick is incorporated in (3.1) by weighting the diagrams 



with quark loops by factor 1/4. The prediction has no free parameters since the pseudoscalar 
masses of various tastes [f7|] and hairpin parameters a 2 5^ A [10] have been determined by MILC. 
In the continuum limit the pseudoscalars of various tastes are degenerate and a 2 8y A — > 0, so the 



lightest intermediate state in ( |3~1| ) is 7TT] . The taste breaking at finite a makes possible the inter- 
mediate state with mass 1M K and is responsible that the effective mass f B SChPT (t) is close to 
2M K , as shown in Fig. ||. The effective mass of B schPT (t) qualitatively agrees with the effective 
mass obtained by MILC [f7|] and UKQCD [Jsj] indicating that the lattice correlators are dominated 
by two-pseudoscalar intermediate states and are significantly affected by the taste breaking. 
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coarse MILC lattice fine MILC lattice 




Figure 4: The squares represent the effective mass of the scalar correlator from MILC staggered simulation 
with 2+1 dynamical quarks. The triangles present the effective mass of B schPT (t) (3.1) at t = 20 for coarse 
MILC lattice and at t = 30 for fine MILC lattice. Mass is plotted in units of r\ ~ 1.6 GeV 1 like in [01 . 



4. Simulations with staggered sea quarks and chiral valence quarks 

The simulations with chiral valence quarks on the available staggered MILC configurations 



[]llj, []_2p present an appealing possibility. However there is no unique recipe how to match valence 
and sea quark masses at finite a and some features of partial quenching always remain. 

The effect of mixed quark actions in simulations with 2 + 1 staggered sea quarks can be at- 
tributed to the bubble contribution fBtl within Mixed ChPT of [fll 1 



B MCHPT (n) _ B 2y J _4 1 l__ (k 2 +Ml Ji )(P + Ml) 

IP) »oL\ 3 {k+p )2+MZ al>val (tf+JK*,,*,) 2 * 2 + ^ + §Mj, (4 " 1} 
1 i 1 i 1 

+ 2 ... ,„ ... + 



The only unknown parameter a 2 AMi x =M 2 al sea —Bo(m sea +m va i), sea = u,s, gives the taste breaking 



in the pion mass composed of one valence and one sea quark Q13|], while the other input parameters 



have been determined in [|7|, |llj]. Fig. |5| shows that the scalar correlator can be negative for light 
u/d quarks if valence and sea quark masses are tuned by matching M va i, V ai = M %i (used by LHP 
(jllj]), while it is positive if M va / )Va / = M %1 are matched j5|]. Comparison of the point-point lattice 



correlators and B MChPT in case of matching M va i va i=M 7t5 offers a possibility to determine a 2 &Mix- 
5. Conclusions 

We determined the mass of the lightest qq scalar meson with 1 = 1 using the simulation with 
two dynamical Domain Wall Fermions. The exponential fit of the scalar correlator is justified in 
this case giving the mass 1.58 ±0.34 GeV. In general the scalar correlator receives also a sizable 
bubble contribution B(t), which is the intermediate state with two pseudoscalar mesons. We pro- 
vide analytical predictions for B(t) for simulations that use (partial) quenching, staggered fermions 



'Our result for B MChPT (t) agrees in the limit t °°, o 2 Amlx —* and m u = = m s with the result of [p^], who 
study mixed quark actions but not the staggered sea. 
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Figure 5: The bubble contribution B MCIlPT (t) (4.1) for simulation with chiral fermions on coarse MILC 
configurations with sea quark masses am u j^ = 0.01 and am s = 0.05. It is plotted for two choices of tuning 
m va [ with m sea and it depends on the value of a 1 k.Mi% , which is varied in the reasonable range [|[] . 

or mixed quark actions. Our predictions for B(t) within relevant versions of ChPT are free of un- 
known parameters, expect for a parameter a 2 &Mix in case of mixed quark actions, which could be 
determined from the scalar correlator. We find that B(t) is sizable at small quark masses and is 
significantly affected by the unphysical approximations that are employed in simulations. It can 
render the negative sign and unphysical effective mass 2M n in the scalar correlator. Our predictions 
for the bubble contribution will be needed in order to extract the scalar meson mass from the scalar 
correlator in the future simulations. 
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